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تحلیل احتمالاتی نشست آنی شالوده‌ی مستطیلی بااستفاده از نرم‌افزارهای مختلف قابلیت اعتماد * 
هادی فتاحی(٩‏ فاطمه جیریایی شراهی( سام مهرنیا ۳ 


چگیفه بررسی نشست در پی‌های سطحی با برعهده داشتن وظیفه‌ی انتقال نیروهای سازه به زمین» یکی از عوامل مهم در طراحی سازه‌ها 
می‌باشد. پارامترهای نحاک به‌دلیل ماهیت غیرهمک نآن. دارای عدم قطعیت‌های زیادی هستند که صرف نظ رکردن ا زآن‌ها در تخمین میزان 
نشست, ایراداتی را در تحلیل پایداری سازه وارد می‌کند. ازاین‌رو, امروزه استفاده از روش‌های قابلیت اعتماد با توانایی پذیرفتن عدم قطعیت‌ها 
و تحلیل پایداری مورد اطمینان مسائل ژلوتکنیکی. مورد استقبال قرا رگرفته‌اند. در این تحفیق به‌منظور تخمین احتمال شکست ناشی از نشست 
آنی یک شالوده‌ی مستطیلی سطحی, از روش‌های مرتبه اول و مرتبه دوم قابلیت اعتماد و هم‌چنین شبیه سازی موئ تکارل و کمک گرفته شد و 
تحلیل‌های مربوط به آن در نرم‌افزارهای 17و 01891 انجام شد. نتایج به‌د ستآمده, حاکی از بالا بودن احتمال شک ست مربوط به د شست 
شالوده کت است. روش‌های مرنبه دوم و شبیه سازی مون تکارلو به‌دلیل حطی نبودن رویه‌ی حالت حدی, نتایج دقیق‌تری را ارائه می‌دهند. 
به‌علاوه تحلیل حساسیت متغیرهای تصادفی براساس بردارهای اهمیت ذشان داد که بار وارد شده بر شالوده نقش مهمی در احتمال شکست 


دارد و پس ا زآن ضریب پواسون متفیرتصادفی بااهمیتی است. 


واژه‌های کلیدی احتمال شکست. نرم‌افزار 387 نرم‌افزار 601856 نشست پی‌های سطحی. 


۱ ۱ 
۱۸۵1۱۵۱01۱۱۲ عبا0ز2۳ ۲۷ 


۵/۵ .5 6 ۲۰ نطم)۲2 ,۲۱ 


ار ۱۰۱۱ 
00 1۶ 10 1۸۱۵ ۶۳۱۵۵۲۵۲۵۶۰ ]0 جع ۱6 ۱ ۵60۲5 11۳0۲۱۵۲ 16 ]0 0۳0۵ 5 ,8۲0۱۵ 1۱۵ 1۵ 0۳65 
6 651111101118 ۱ ۱۷۵/۱۲۵1 0۴ 01۱۱۱59/0۲ 11۱6 رک0۱6۵۱۵ ۵ ۱۵۷ 00۱۵25 50 ,1۵۵/۵ 01۱5010۳۵۱6 01 
0۴ ۱56 1۱۱6 ,1۱۵۳۵/۵۲۵ ۵اه ماک ۵ کاع06۵ ۱6 ۳ 61۱025 0۳۵۵65 ,0۷۱۵۵۵1 0۴ 561116111۵۲1 
۵ 0 0۵5۱ وهای ع۲۵]۱۵ ۵۱۱4 ۷۱6۵1۵6۶ 066601 10 ]۵91 1۱6 ۷۲/۰ ۱۱6۵0۵5 و۲66 
0۴ هی 10 0۵ 0۵۵۵ چاااراهرا2۲۵ ۱6 6513۱۱۵۱۵ 10 0۲6۲ 1 ,اک 15 ۱ ۱۵۵6 ۱۷۷۵۱6۵۱۱۵۵ 0۲۵ 1951/65 
۵ ۱۷/0۲۱/۵۶ که ]۱۷۵ 6 ۱6۳۵85 طز۲۵]۱۵۳ 66600-0۲۲ ۵۳00 ۲5-0۲۵۲ ص0۵0 ۲۵616۵ 
6 1116 .۹0۵1۱۵۵۲۵ /۵6) ۵0 ۳ 1 ۵/۵۳۱۵0 ۱۵۵۲۵ ک0۳065 ۲۵/۵۷۵ 1۱6 0۱ ۱6۵۵ ۷۷۵۳۵ 5111۵۵015 
015 56601۱۵-07۵۲ 1116 9.5۵(۰) 661116111۱611۴ 1۱60۱۵۵101 10 ۲6۵۱۵۵ 0۳۳۵۵۱۵۲۵ ]2۳006۵ حلع 6 11۱0166016 
6 ۱۱۱ 1/۱6 0۴ ۱۱۵۱-۱۵۵۲1۵۲ 16 10 6 ۲۵۵۱۵6 066۵۳۵۵ 1۱۵0۳۵ 2۲۵۷۵ 53۵005 ۵۲0 ۷۹۵۱۱۱۵ 070 
۵ ۷6۵0۲۶ ای 01 هر عع۷۵۳۱۵۵ ۲۵/۱۵۵۲ 0 امه 6۵۲ ۱6 رص00 ۱ ۲۵۵۵۵۲ 
0 ۶ 015901 ۱۱۵ ۱۱۵۴۲ 00 ,060۲6 زا۳۳۵۵ 6 ۲ ۲۵6 1۱۳۴0۳۱۵۲۱۶ 6۱ 66 ,(0) ککت۳ک 601۱1661 1116 17101 

9 1 


۵۰ 1201088 ص10 ممه یاو رمته)۵1و دنه رمتهسگمو ۴1 ,عبانم گه نازمد۲۸۵۵ ۲۷۷۵۸۲۵ مک 


* تاریخ دریافت مقاله ۹۹/۵/۲۲ و تاریخ پذیرش آن ۹۹/۷/۲۷ می‌باشد. 
(۱) نویسندة مسئول, دانشیار» دانشکده‌ی مهندسی علوم زمین. دانشگاه صنعتی اراک. طع.بعلمیع ۵ تطم))ه1.ظ خاتمصصرز 
(۲) دانش آموخته دانشکده‌ی مهندسی علوم زمین, دانشگاه صنعتی اراک. 


(۳) دانش آموخته, دانشکده‌ی مهندسی علوم زمین, دانشگاه صنعتی اراک. 


۱۸ تحلیل احتمالاتی نشس تآنی شالوده‌ی مستطیلی باستفاده از نرم‌افزارهای... 


مقدمه 
تمامی سازه‌هایی که برروی زمین بنا می شوند از جمله 
سساشنمان فان پل‌ها و اک بوهاه از دو نکن سسازدق 
فوقانی و سازه‌ی زیرین تشکیل می شوند. بخش سازه‌ی 
زیرین که به‌عنوان پی یا فونداسیون از آن یاد می‌شسود» 
به‌عنوان حائل بین سازه‌ی فوقانی و زمین تکیه‌گاه عمل 
می‌کند. به‌نحوی که تنش‌های بیش از حد و نیز 
نشست‌های اضافی ایجاد نشود. در طی فرآیند 
شکل‌گیری نهشته‌های رسوبی قبل از هرگونه بارگذاری» 
ناهم‌سانی ذاتی در خاک‌ها ایجاد می‌شود هم‌چنین در ثر 
اعمال سربار به‌علت کرنش‌های ایجاد شده‌ی نا شی از 
تغییرات تنش, تغییر جهتی در ذرات خاک ایجاد می‌ شود. 
که این امر ناهم‌سانی القایی را ایجاد می‌کند. بیشترین 
تأثیر این ناهم‌سانی برروی پارامترهای الاستیسیته‌ی 
حاک است. یکی از مهم‌ترین کاربردهای این پارامترهای 
الا ستء سیته‌ی خاک تخمین نذشست و برر سی ظرفیت 
باوبری پی‌هاق ستطحی اسست[1] عوامل .ذکرشسدهق 
برخی عوامل دیگر باعث وجود عدم قطعیت در 
پارامترهای خاک می‌شود که محاسبات و تحلیل‌های 
قطعی را دچار خطا می‌کند [2]. به‌عبارتی خاک و سنگ 
از جمله محیط‌هایی هستند که برخی ویژگی‌های آنها در 
هر دو نقطه از آن متفاوت است و باعث ایجاد ضریب 
تغییرات می‌شود. دراین صورت چاره‌ای جز یافتن راهی 
برای ورود عدم قطعیت‌ها نیسست. به‌طوری که چشسم 
پوشی از آنهء نتایج تحلیل‌ها را از واقیت دور می‌کند و 
صمیم‌گیری در ستی برا ساس آنها انجام نمی شود. در 
سال‌های احیر جهشی برای استفاده از نظریه‌های قابلیت 
اعدماد برای مدل‌کردن عدم قطعءیت ها در طرح های 
مهند سی به‌وجود آمده‌ا ست [3,4]. در طراحی برا ساس 
قابلیت اعتماد به‌جای استفاده از بردار مقادیر پارامترها 
به‌صورت قطعی " [م,...,2,:] < ۵6 توزیع آماری 
آنها درنظر گرفته می‌شود. به‌عبارتی درصورت مشخص 
بودن توابع چ گالی احمال برای پارامتر های بار و 
تلاوت من فان از طرالفی پر ایستاش فاشت: ایا 
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استفاده نمود. محاسبه‌ی احتمال شکست برگرفته شده از 
تحلیل قابلیت اعتماد به‌کمک تابع حالت حدی انجام 
می‌شسود. حد تعادلی اين تابع (0 < (006) مرز بین 
منطقه‌ی ایمن و شکست است [3]. شکل (۱) تابع 
چچگالی استمال (6)10 و چگونگی تفکیک قسمتی از 
سطح زیر نمودار را به‌عنوان احتمال شکست نشان 


می‌دهد. 


662( 


شکل ۱ نمایی از تعیین احتمال شکست بر توزیع احتمالاتی تابع 
حالت حدی [6] 


پیچید گی محاسبات احتمال شکست در روش‌های 
تحلیلی که از تعبین مساحت زیر نمودار و انتگرال بسیار 
پیچیده منتج می‌شود. سبب شده‌است که در برحی 
روش ها | ندیس قابایت اعتماد (8) جایگزین این 
محاسبات پیچیده شود [5]. 
)۱( ( )۵ ۲۳۵۲ 
به‌طوری که ۵ تابع توزیع تجمعی | ستاندارد نرمال 
افو ععاشیی این آلدیس» ها از تفای ال 
هرد به فشسای استفاوداری پمال اقال پیداامی کنند: 
اندیس قابلیت اعتماد در روش تحلیلی از طریق 
محاسبه‌ی ممان‌های مرکزی تابع حالت حدی به‌دست 
فیس ایا اندیسی که توسط هاسغفر و لیند (-۲۱۵90]6 


سال سی و سوم شمارة سهء ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی- سام مهرنیا 


9 ( - 70( - 06 /جنود < 8 

که در آن 26 بردار متغیرهای تصادفی (6» ۲ بردار 
میانگین متغیر هاء » ماتریس کوار یانس و "3 دامنه‌ی 
شکست است. در اینجا 8 از طریق یک مسئله‌ی 
بهینه سازی به‌د ست خواهد آمد. روش ممان‌های آماری 
ماذند روش های مردبه‌اول و مردبه‌دوم قابایت اعتماد 
ازجمله روش‌هایی هستند که بااستفاده از اندیس قابلیت 
اعتماد به تخمین احتمال شکست می‌پردازند. 

م طااعاتی که درمورد کاربرد قابایت اعتماد در 
برحی مسائل ژئوتکنیکی. در گذشته توسط محققان 
انجام شده‌است. به شرح زیر است: 
بابو و هم کارانش تحایل قابلایت اعتمادی برمب نای 
محاسبه‌ی فشار مجاز فوندا سیون‌های سطحی متکی بر 
خاک‌های رسی سخت انجام دادند [7]. 
گریفث و فتون در سال ۲۰۰۷ با تلفیق تحلیل احتمالاتی 


و روش اجزای م‌ حدود. تراوش آب در زیر سد ها 


ظرفیت باربری و ند ت پی‌ها و پایداری پایه‌های 
معدنی را مورد بررسی قرار دادند [8]. 

تفر کت که بایان نت ها ان اما رد 
مهندسی ژئوتکنیک هستند مطالعات متعددی برروی 
موارد عدم اطمینان در ظرفیت باربری پی‌ها در خاک‌های 
دانه‌ای و چ سبنده انجام داده‌اند. نظریه‌ی اجزای محدود 
تصادفی از جدیدترین تئوری‌های آنها در سال ۲۰۰۸ 
می‌باشد که از ترکیب روش اجزای محدود و روش 
شوک قازان خاضا فتوایست: [9] 

لو و همکاران . تحلیل قابلیت اعتماد را برروی شکست 
گوه‌ای تو سعه دادند و ارتفاع گوه را تغییر دادند و چهار 
حالت شکست گوه‌ای را بررسی کردند [10]. 

سانجی و کیچی در سال ۲۰۰۹ آنالیز قابلیت اعتمادی 
انیجام داد ند و روان‌گرایی خاک را براساس داده‌های 
آزمایش نفوذ استاندارد مورد بررسی قرار دادند [11]. 
آقای ناگ و همکاران در سال ۲۰۱۱ پی‌های سطحی در 
خاک چ سبنده را در شرایط مختلف بار گذاری و 


شبیه سازی کردند. به‌طوری‌که زاویه‌ی | صطکاک خاک 


سال سی و سوم. شمارة سهء ۱۳۹۹ 


بارهای افقی و عمودی را به‌عنوان متغیرهایی با توزیع 
نر مال در مطااعه در نظر گرفتند و تابع عملکرد را با 
معادله‌ی مایرهوف به‌وجود آوردند. نتایج ذ شان داد که 
احتمال ‏ شک ست در پی‌های سطحی متأثر از تنوع بار 
عمودی. زاویه‌ی اصطکاک و بار افقی می‌باشد [12]. 
چینگ و همکاران در سال ۲۰۱۱ مطالعه‌ی آنالیز قابلیت 
اعتماد ظرفیت باربری پی‌ها برروی لو دومن تسیک را 
بااستفاده از مطالعه‌ی موردی منطقه‌ای در تایوان انجام 
دادند و به این نتیجه ر سیدند که شاخحص مقاومت زمین 
شنگ» دی پارامتر تاثیر گذار هستند. [13], 

چین لونگ و بک. قابلیت اعتماد شمم‌ها یا المان‌های 
انتقال بارهای سطحی را به عمق زمین بااستفاده از 
احتمالات تحلیل کردند [14]. 

لوانیس و فیایپ با به کارگیری روش مو نت کارلو, 
پایداری داتخلن دیوار خاکی تقویت‌شده را بررسی کردند 
[15]. 

جوهری و سبزی. بااستفاده از روش‌های مختلف عددی 
موفق به پیش‌بینی میزان ذدشست شالوده شدند و آنالیز 
قابلیت اعتماد میزان نشست را نیز انجام دادند [16]. 
نظرزاده و سسربیشه تحلیل قابلیت اعتماد شالوده‌ی 
سطحی را انجام دادند و بر این اساس برای متغیرهای 
ورودی توابع توزیع احتمالاتی را درنظر گرفتند. آنها در 
این آنالیز از روش شبیه‌سازی مونت کارلو استفاده کردند 
[17]. 

شکیر » ظرفیت باربری یک شالوده‌ی مربعی را محاسبه 
کرد و با ستفاده از روش مونت کارلو تابع حالت حدی 
آن را مدل کرد و به تحلیل احتمالاتی آن پرداحت [18]. 
در این پژوهش سعی شد بااستفاده از توابع توزیع 
متغیرهای تصادفی در مسئله‌ی ند ت شالوده‌ی 
مستطیلی سطحی و هم‌چنین روابط تعریف‌شده یک 
آ نالیز قابا یت اعت ماد ان جام شود. به‌این‌منظور. از 
روش های مختلف قابلیت اعتماد مانند روش قابلیت 


اعد ماد مرد به‌اول» روش قابایت اعد ماد مرذ به‌دوم و 
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شبه سازی مونت‌کارلو برای تخمین احتمال شکست و 
این تحایل ها دوه نرم‌افزار های 1 و ۵165 که ابزار 
قدرتمندی برای انجام تحلیل های احتمالاتی هستند 


صورت گرفت. 


و هی یت هن من ی بای 
روش‌های ارزیابی قابلیت اعتماد که در نهایت به تخمین 
احتمال شکست سیستم منجر می شوند. به‌طور تقریبی 
نسبت به روش‌های تحلیلی. هدف خود را دنبال می‌کنند 
تا به قضاوت درمورد پایداری سیستم‌ها بپردازند. روش 
مرتبه‌اول قابلیت اعتماده روش مرتبه‌دوم قابلیت اعتماد و 
شبیه‌سازی مونت‌کارلو از جمله‌ی اين روش‌ها هستند که 


به آنها پرداخته خواهد شد. 


روش مرتبه‌اول قابلیت اعتماد 
تابع حالت حدی در مقدار صفر آن (0 < (0:000) که 
رویه‌ی حالت حدی نام دارده ممکن ا ست به شکل یک 
خط راست تا یک منحنی با درجه‌ی بالا باشد. در روش 
ممان‌های آماری» تخمین حجم احتمال مشسترک 
قرارگرفته در منطقه‌ی شکست. به‌منظور محاسبه‌ی 
احتمال شکست. با فرض‌هایی ذسبت به شکل رویه‌ی 
حالت حدی انجام می‌شود. در روش مرتبه‌اول. که 
ترکیبی از روش‌های تحلیلی و تقریبی است. این فرض 
مبنی بر خطی بودن رویه‌ی حالت حدی می‌باشد و 
بااستفاده از بسط درجه‌اول تیلور خطی سازی می شود تا 


منجر به سادگی محاسبات پیش‌رو شود [4] 


600 2 60( + ۷۵0۵۵0۵6-۴ (۳) 


در این رابطه ۶ بردار متغیرهای تصادفی 
1 
نقطه بر رویه‌ی حالت حدی است. بنابراین پس از انتقال 
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تحلیل احتمالاتی نشست آنی شالوده ستطیلی بااستفاده از نرم‌افزارهای... 


استاندارد نرمال و محاسبه‌ی شاخص قابلیت اعتماد. 
احتمال شکست با رابطه‌ی (۱) به‌دست خواهد آمد [4]. 


روش مرتبه‌دوم قابلیت اعتماد 
در اکثر مواقع تابع عملکردی که به‌عنوان یک متغیر 
تصادفی وابسته از متغیرهای تصادفی مستقل تأثیر 
می‌پذیرد. از یک خط را ست پیروی نمی‌کند و متنا سب 
با میزان پیچیدگی آن, مطابق با یک منحنی در جه‌دو یا 
بالاتر ا ست. در روش مرتبه‌دوم قابلیت اعتماد از رویه‌ی 
حالت حدی حول محدمل‌ترین نقطه‌ی آن تقر یب 
درجه‌دوم زده می‌شسود. به‌عبارتی؛ به‌جای تقریب زدن 
رویه‌ی حالت حدی با یک صفحه يا ابرصفحه با یک 
سهمی يا ابرسهمی جایگزین می‌شود [4]. 
ی 
08 
۳ج ( ۷ ب 2 ور ,ژارتم 
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0428 
021 021 


یروا بر 
2 
که در آن :26 متغیرهای تصادفی و :5 برآمدهای آنها 
هستند و هم‌چنین 5 بیانگر محتمل‌ترین نقطه اسست. 
بریتنگ (ع0اع:8) رابطه‌ی زیر را برای محاس‌به‌ی 
احتمال شکست با تقریب مرتبه‌دوم ارائه داد [4]: 
(0) 1 
که در آن :۷ نشان‌دهنده‌ی انحنای اصلی تابع حالت 
حدی در نقطه‌ی طراحی است به صورتی‌که اگر تابع 0 
متعین دنه با شتلم: (ابض) انا دارد بات این تثاز است 
که این انحناها محاسبه شوند [4]. 


روش شبیه‌سازی مونت کارلو 
روش شبیه‌سازی مونت کارلو ازجمله روش‌های عددی 
است که در آن هر متغیری به‌طور هم‌زمان طبق قوانین 
تعیین‌شده. تغییر می‌کند. ساده‌ترین قانون این است که 
برای یک تحلیل پایداری داده‌شده هر متغیر تنها یک 
مقدار انتخاب‌شده‌ی تصادفی از توزیع اندازه گیری 
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دی که مهار اون مان مها ی باق تاه 
داش شین رفاووضل سس که وو مدا ماخ 
حدی نمود پیدا می‌کند به|زای این مقادیر انتخاب‌شده 
محاسبه می‌شود و به‌این ترتیب یک آزمون از رفتار سیستم 
انجام می‌شود. مونت‌کارلو این آزمون‌ها را به تعداد بسیار 
زیاد تکرار م‌کندبه‌طوری‌که توزیع احتمالاتی تایع سالت 
حدی مشخحص شود. به این ترتیب احتمال شکست 


به‌شکل زير محاسبه می‌شود [19]. 
)0( 


(8660 
آّ 


تن 
0 


که در آن (0 > 166 تعداد آزمون‌هایی است که در 
آنها 0 > 6 است. تعداد آزمون‌ها آا متناسب با دقت 
احتمال شکست مطلوب می‌باشد. از مزایای مهم این 
روش عدم وجود هیچ محدودیتی در نوع و شکل توزیع 
احتمالاتی متغیرها و نیز تابع حالت حدی می‌باشد [19]. 


بردارهای حساسیت 

متغیرهای تصادفی. هر کدام به نوعی بر پایداری سیستم 
و احتمال شکست آن تأثیر می گذارند. آگاهی از میزان 
تأثیر گذاری متغیرها و علاوه‌بر آن یافتن با اهمیت‌ترین 
آنها؛ کمک زیادی به طراح و تدارک اقدامات کنترلی 
برای پایین آوردن احتمال شکست می‌کند. باا ستفاده از 
ماهیت تغییرپذیری متغیرها می‌توان تأثیر آنها بر شاخص 
قابایت اعتماد و احتمال شسکست را مطااعه کرد. 
بردارهای حساسیت آلفاء بتاء دلتاء اتا و کایا از جمله 
بردارهای حساسیتی هستند که با ارائه‌ی اهمیت متغیرها 
در قالب یک بردار به مقایسه‌ی متغیرها نسبت‌به یکدیگر 
می پرداز ند. بردار حساسیت آلفا؛ اهمیت هریک از 
متغیرها را براساس سهمی که در محاسبه‌ی واریانس 
حدی از این جهت که در محاسبه‌ی احتمال شکست 
دخیل است بیان‌کننده‌ی این اهمیت است [14]. 


2 رد ...دوه + ۵۶2) ۱۳۵۱ - [۷]6 
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حدی به اندازه‌ی 03 یعنی بردار یکه‌ی آن است. هرچه 
]تشن تاشست مانگر آن اسنتت که اهست وبا نو 
ات بعااشتی وه نیز در احمال شکست: بات کذار اشت. 
مقدار مثبت » نشان می‌دهد با افزايش آن» متغیر ۲ 
افزایش پیدا می‌کند و مقدار منفی بیانگر آن است که با 
افزایش آن متغیر :۳ کاهش می‌یابد. در بردار حساسیت 
دلتا تغییرات شاخحص قابلیت اعتماد نسبت‌به میانگین 
متغیرها سنجیده می‌شود. این نشان‌دهنده‌ی اهمیت 
میانگین متغیرهای تصادفی است [4ا1. 


8 ۳ 0۷ ٩ (۸) 


هم جنین در بردار اهمیت اند تغییرات شاخحص 
قانایت اعتمان شنت به بر اکنط کی تیه پروستی 
می‌شود. 

0 
)0 ره < را 
نشست آنی شالوده‌ی سطحی 

دو نوع نشست اصلی می‌تواند برای شالوده‌های سطحی 
رخ دهد؛ نخست نشست آنی يا الاستیک و دوم نشست 
تحکیم. نشست آنی شالوده بلافاصله بعد از احداث 
شالوده رخ می‌دهد. ذشست تحکیم یک ذشست تابع 
زمان می‌باشتد و مطابق شوری تحکیم علت وقوع آن 
خروج آب حفره‌ای از میان حفرات خاک به‌علت 
اضافه‌فشار تحمیل‌شده بر آنها می‌باشد. نشست کلی یک 
شالوده مجموع ذش ست‌های ۳ (الا ستیک) و تحکیم 
می‌باشد [20]. از آنجا که در محاسبات مربوط به نشست 


آنی از پارامترهای الاستیک خاک استفاده می شود به آن 
نشست الاستیک نیز گفته می شود. محاسبات ذشست 
آنی يا کو تاه مدت در مورد خاک های دا نه‌ای و نیز 
خاک‌های چسبنده با درجه‌ی اشباع کمتر از ۸۰ در صد 
کارپرد دارد [21]. شکل (۲) ترسیمی از شالوده‌ی 
مستطیلی سطحی و متغیرهای تأثیر گذار در میزان نشست 


را نشان می‌دهد. 
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شکل ۲ ترسیمی از شالوده‌ی مستطیلی سطحی و متغیرهای 


تأثیرگذار در میزان نشست 


معرفی اجمالی نرم‌افزارهای به کار گرفته‌شده در 
این تحقیق 
نرمافزارهای متعددی قابلیت انجام تحلیل‌های احتمالاتی 
و ارزیابی ریسک و هم‌چنین آنالیه سا پنسا قار نز 
در این تحقیق از دو نوع آنها برای انجام تحلیل‌ها کمک 
گرفته شده‌است. در ادامه به تعریف دو نرم‌افزار 1*7 و 


۵۳8 پرداخته خواهد شد. 


نرم‌افزار 11 

قابلیت اعتماد و ریسک را دارد. پارامترهای سیستم 
ازجمله متغیرهای تصادفی و قطعی. ورودی‌های این 
برنامه هستند. در این برنامه امکان استفاده از تعداد زیادی 
متغیر ورودی و ارزیابی حساسیت و اهمیت آنها نسبت‌به 
یکدیگر وجود دارد و در آن می‌توان از توزیع‌های 
مختلفی مانند نرمال. گامبل بتاه نمایی لاگ نرمال و رایلی 
برای تعریف رفتار متغیرها استفاده کرد. انجام شبه‌سازی 
مونت کارلو. روش مرتبه‌اول قابلیت اعتماد و روش 
مرتبه‌دوم قابلیت اعتماد از دیگر قابلیت‌های برنامه است. 
با انجام اين آنالیزهاء مقداری از متغیرها که به‌ازای آنها 
هم‌چنین این برنامه با برنامه‌هایی مانند ۸۵۳ 
5 دنمدطض وعع۵28ع۵ر تلض]1 ۷۲۸ می‌تواند 
لینک شود. این نرم‌افزار را می‌توان به‌صورت رایگان از 
طریق لینک (107191006.02/50117۷276/71) دریافت کرد. 


نشریهةُ مهندسی عمران فردوسی 
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نرم‌افزار ن08) 
نرم‌افزار ریسک یک سیستم نرم‌افزاری قوی درزمینه‌ی 
تحلیل وضعیت‌های کاری و فنی است که تحت تأثیر 
ریسک قرار دارند. اين نرم‌افزار تکنیک‌های خود را در 
قالب بسته‌ی نرم‌افزاری صفحه‌گسترده‌ی استاندارد 
نایکروسافت اکما اوانه موه و ماصرست نی اون 
۳ برروی اکسل نصب می‌شود. سیستم نرم‌افزار ریسک 
تمام ابزارهای مورد نیاز برای تعریف مجموعه, اجرا و 
نمایش نتایج تحلیل ریسک را ارائه می‌دهد و در وضعیت 
مشابه‌با مت هو عملگرهای اکسل عمل عی کند. بااستظاده 
از نرم‌افزار ریسک و اکسل هر وضعیت ریسکی از 
تجارت و کسب و کار گرفته تا علوم و فنون مهندسی 
می‌تواند مدل شود. توزیع‌های احتمالی ارائه‌شده در آن 
امکان طبقه‌بندی نمودن تقریبا تمامی انواع عدم قطعیت‌ها 
رخا سار ای مت ترجه سجن یس 
بااستفاده از روش‌های شبیه‌سازی مونت کارلو و یا ما1 
۵6 و تولید اعداد تصادفی. مدل را بارها اجرا 
کرده و پارامترهای پروژه را تحلیل می‌کند و درنهایت 
مقادیر خروجی‌های مشخص‌شده و توزیع احتمال آنها را 
ارائه می‌کند. این نرم‌افزار به‌منظور انجام سریع‌ترین 
شبیه‌سازی ممکن و بااستفاده از تکنیک‌های نمونه‌گیری 
رف طراننی بات خی اهر ال که شیباا با 
وضعیتی مواجه شدید که می‌بایست یک تصمیم و یا 
تحلیل همراه با ریسک انجام دهید. نرم‌افزار ریسک 
می‌تواند به‌منظور بهبود چشم‌انداز تصوير وقایعی که در 
آینده رخ خواهند داد استفاده شود. پس از ورود متغیرها 
و تعریف توزیع احتمالاتی آنها پرنامه آماده است تا با 
تولید اعداد تصادفی و بااستفاده از روش‌های احتمالاتی 
خود مدلی بسازد و برای متغیر خروجی یک توزیع 
احتمالاتی ارائه دهد. درنتیجه تحلیل‌ها براساس رفتار 
توزیع ارائه‌شده انجام ین کف از مزیت‌های این نرم‌افزار 
تولید شکل‌های گرافیکی پیرامون ورودی‌ها و نتایج است 
که باعث فهم بهتر و سهولت در تحلیل‌ها می‌شود. این 
نرم‌افزار قابلیت کاربرد در هر زمینه‌ای را که عدم قطعیت 


ذر آن تعابان است دار 
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تحلیل و نتایج 
در ادامه به بررسی ند ِ آنی یک شسالوده‌ی 
انعطاف‌پذیر مستطیلی به‌منظور برآورد ایمنی آن در یک 
م سئله‌ی ژئوتکنیکی پرداخته خواهد شد. در تحلیل‌های 
انجام‌شده روش‌های مرتبه‌اول و دوم قابلیت اعتماد و 


شبیه‌سازی مونت‌کارلو به‌کار گرفته شده‌است که در 
نرم‌افزارهای 1۲ و ۵185 مدل‌سازی شده‌اند. 

مسئله‌ی نشست شالوده‌ی مستطیلی باتوجه به 
رابطه‌ی زیر شامل متغیرهای تأثیرگذاری است [22] که 
جدول (۱) به شرح پارامترهای قطعی این مسئله 


می‌پردازد. 
ِِ 2 
۰ 3 ( را 2 ,)سس ووظ0.5 - ۸1۷ 


1-۷ 


برخی از اين متغیر ها به‌دلیل وجود عدم قطعیت 
بالای ویژگی های ژئوتکنیکی. به‌عنوان متغیرتصادفی 
درنظر گرفته‌نده‌اند. جدول (۲) بیانگر پارامتر های 
احتمالاتی متغیرهای تصادفی می‌باشد. به‌علاوه شکل‌های 


۳۳ 


(۳) تا (۵) توزیم احتمالاتی متغیرهای تصادفی را نمایش 
وارده. ضریب پواسون و مدول الاستیسیته به‌عنوان 


متغیرهای تصادفی | علامت اختصاری 
بار وارده («ک) وه 
ضریب پواسون ۷ 
مدول الاستییته (۷۲۳۵) ,۴ 


شده‌اند. 
جدول ۱ پارامترهای قطعی و مقادیر آنها 
پارامترهای قطعی علامت احتصاری | مقدار پارامتر 
عرض شالوده (3)۳ ۳9 
طول شالوده ۳ ۰ 
عمق جایگزینی شالوده (1)۳ ۳ 
ضخامت طقه ()۲۱ ۲۰ 
تعداد گوشه‌های شالوده 1۳ 3 
ضریب کف آ ۷۳ 
ضریب تأثیر ۳ ۰/۱۸۹ 
ضریب تأثیر 1 ۹0۵ 
توزیع احتمالاتی | میانگین | انحراف معیار 
نرمال ۲۳۸۰ ۶۰ 
نرمال ۳۵ ۰/۸ 
نرمال 0۰ ۳/۵ 


1 
۰ 
3 
4 


شکل ۳ توزیع چگالی احتمال متغیر بار وارده 
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شکل ۶ توزیع چگالی احتمال متغیر ضریب پواسون 


تابع چگالی احتمال 


شکل ۵ توزیع چگالی احتمال متغیر مدول الاستیسیته 


روش‌های تحلیل قابلیت اعتمادی چون روش‌های 
مرتبه‌اول و مرتبه‌دوم مقدم بر تخمین شاخص قابلیت 
اعتماد و احتمال شکست. نیاز به یک تابع حالت حدی 
به‌منظور ارزیابی عملکرد تابع اصلی در دامنه‌ی متغیرهای 
تصادفی دارند. یک رویه‌ی مشخص برای ساخت توابع 
حالت حدی از توابع صریح, ترکیب آنها با حد آستانه‌ی 
متغیر م ستقل می‌با شد. انتخاب حد آ ستانه به‌عنوان یک 


نشریه مهندسی عمران فردوسی 


عخالتن ادابه مسنسیر فانایت اختماف را تحت تاتیر قرار 
می‌دهد. میزان ذشست شالوده به‌عنوان متغیر مستقل در 
رابطه‌ی (۱۰) با ترکیب حد آستانه‌ی نشست (۵۰ 
)۱۱ +نهنر(۵1) < (6)») 

روش های مرتبه‌اول و مرئبه‌دوم قابلیت اعتماد در 


سال سی و سوم, شمارهُ سه, ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی- سام مهرنیا 


تخل را به‌منظور یافتن نقطه‌ی طراحی به‌کار 
گرفته‌اسست. باتوجه به شکل (۱) در گام سوم از 
جستجوی این الگوریتم نقطه‌ی طراحی ارائه شده‌است. 
در این شکل نمودار آبی‌رنگ فاصله‌ی نقطه‌ی مورد تا 
نقطه‌ی طراحی و نمودار قرمزرنگ. فاصله‌ی نقطه‌ی 
مورد نظر تا مبداً را نشان می‌دهد. درنتیجه فاصله‌ی 
نقطه‌ی برآوردشده در گام سوم تا مبدأء نشان‌دهنده‌ی 
شاخخص قابلیت اعتماد می‌باشد. 

خروجی محاسبات انجام‌شده و احتمالات در هر 
یک از روش‌ها در جدول (۳) آورده شده‌است. احتمالات 
تخمین زده‌شده پیرامون حالت شکست در روش‌های 
مرتبه‌دوم قابلیت اعتماد و شبیه سازی مونت‌کارلو بسیار 
به‌هم نزدیک هستند. درحالی که روش مرتبه‌اول قابلیت 
افقنای مفلار انشمال. شکست بالا فرش را سیخ به سایر 


روش‌ها تخمین زدها تا به این دلیل که نتیجه‌ی روش 


۳۵ 


مرتبه‌دوم به شبیه سازی مونت کارلو به‌عنوان یک روش 
با دقت ب سیار بالا نزدیک‌تر ااست. می‌توان نتیجه گرفت 
رویه‌ی حالت حدی (0 < (06060) در این مطااعه از 
بنابراین در مورد این مسئله روش مرتبه‌اول برای تخمین 
احتمال شکست با تقریب خطی دجار خطا می‌شود. 
دقت بالا» تعداد بسیار زیادی داده‌ی تصادفی تولید 
می‌کند. هرچه تعداد تکرار شبیه‌سازی‌ها بیشتر شود 
احتمال تخمین زده‌شده به مقدار واقعی آن نزدیک تر 
اشحت: بااستفاده از شکل (۷) که روند هم‌گرایی روش 
شبیه‌سازی مونت کارلو به جواب و کاهش تغییرات در 
۰ داده‌ی تصادفی از متغیرهای تصادفی روش 
مونت کارلو به احتمال ۹/۵۶ هم‌گرا می‌شود. 


۱/۳۵ 


9 


۶ ۵ ۴ ۳ ۲ ۱ 
۱۵۲ مع5 


شکل 1 نمودار هم‌گرایی فاصله تا نقطه‌ی طراحی 


جدول ۳ احتمال شکست محاسبه‌شده با روش‌های قابلیت اعتماد 


‌ ۳ 45 ۱/۵ 
قابلیت اعتما 

تحلیل قابلٍ د‌ ۱9 52007 ۳ ودنه) 

احتمال شکست (/) ۳ ۹/5۹ ۹/۵۶ ۹۹ 

شاخس قابلیت اعتماد ۱/۳۳۷ ۱/۳۰۵ ۱/۳۸ ۳ 


سال سی و سوم شمارة سهء ۱۳۹۹ نشریه مهندسی عمران فردوسی 


ک 


تحلیل احتمالاتی نشست آنی شالوده صستطیلی بااستفاده از نرم‌افزارهای... 


کواریانس احتمال شکست 

1 
08 
06 
04 
02 

0 ۳ 

20 40,0 200 72-0 1,0 0 
تعداد نمونه ها 


شکل ۷ نمودار کواریانس احعمال شکست تسب به افزایش تعداد نمونه‌ها 


هم‌چنین شبیه‌سازی مونت‌کارلو با نرم‌افزار 6۵1551 
خروجی مشابهی را ارائه می‌دهد. در شسکل (۸) تابع 
توزیع فراوانی (026 رسم شده‌است. با محاسبه‌ی 
قسمتی از مساحت زیر نمودار که در سمت منفی نمودار 
قران گر فتاه استمال شسسکشتک:۱ ۰۹/۳ به‌دست 
ال واینیت: 

گاهی توابع حالت حدی دارای روابط بسیار 
پیچیده‌ای هستند که انجام شبیه‌سازی مونت‌کارلو برای 
دست یافتن به دقتی بالا زمان و هزینه‌ی محاسباتی بالایی 
را صرف می‌کند. درحالی‌که روش‌های مرتبه‌اول و 
مرتبه‌دوم قابلیت اعتماد با صرف زمان کمتر به جواب 
مشابهی می‌رسد. در تحقیق حاضر روش مرتبه‌دوم 
نسبت‌به روش مرتبه‌اول خروجی دقیق‌تری را از احتمال 
شکست ارائه می‌دهد و نسبت به شبیه‌سازی مونت‌کارلو 
هزینه‌ی محاسباتی کمتری دارد. 

تحلیل حساسیت متغیرها و تأثیرگذاری آنها با 
جداول و نمودارهای مربوط, دید مناسبی نسبت‌به اهمیت 
آنها در برآورد ایمنی سازه‌ها؛ به طراح می‌بخشد. نتایج 
آنالیز حساسیت برمبنای قابلیت اعتماد باتوجه به شکل 
)٩(‏ بیانگر حساسیت بالای شاخحص قابلیت اعتماد 
نسبت‌به متغیر بار واردشده بر شالوده‌ی سطحی می‌باشد. 
واضح است باری که درنهایت بر شالوده وارد می‌شود در 
میزان نشست آن تأثیر زیادی دارد. اما کم‌کردن این بار 
برای بهبود اوضاع شالوده در مواردی ناممکن است. 
چیزی که اهمیت دارد بهبود بخشیدن پارامترهای محیطی 
است که متحمل بار می‌شود. در این مورد ضریب پواسون 


از اهمیت بالاتری نسبت به مدول الاستیسیته برخوردار 


نشریهةُ مهندسی عمران فردوسی 


انسشت الیقه تتیتجه گر دگراشله راشای شبه سردا اهشیت 
آلفا ات و علتا ارائه شنه‌است: 
شکل (۱۰) نشان‌دهنده‌ی اهمیت متغیرها براساس 
بردارهای اهمیت بثا و کایا است که اولویت تأثیر گذاری 
متغیر ها را نسبت به سه بردار اهمیت دیگر تغییر 
داده‌است. 

جدول (4) به‌طور دقیق‌تر بیانگر اهمیت متغیرها 
براسامن. .هر فا شمیت ۰ ایک ماخ اخست 
محاسبه‌شده ممکن است مثبت پا منفی باشند. این علامت 
در بردارهای مختلف مفهوم متفاوتی از اهمیت متغیرهای 
مورد بررسی دارد. بردار اهمیت آلفا و کاپا با درنظرگرفتن 
علامت منفی برای حساسیت متغیرهای مدول الاستیسیته 
و ضریب پواسون, نشان دادند که اين متغیرها از نوع متغیر 
مقاومت هستند و با افزايش مقادیر آنها احتمال شکست 
کاهش می‌یابد. به‌علاوه با درنظرگرفتن علامت مثبت برای 
متغیر بار وارده. ثابت کردند که اين متغیر از نوع متغیر بار 
است و با افزايش مقدار آن احتمال شکست افزایش 
می‌یابد. باتوجه به اين که بردارهای اهمیت بتا و دلتا عکس 
این رویه را دارند. علامت‌ها به‌طور برآوردشده است. اما 
بردار اهمیت اتا با در نظر گرفتن پراکندگی هر متغیره 
حساسیت متغیرها را منفی فرض می‌کند. به این دلیل که 
با افزايش انحراف معیار متغیرها مقدار احتمال شکت 
افزایش می‌یابد. 
شکل (۱۱) تحلیل حساسیت در نرم‌افزار کلونگه را 
نشان می‌دهد که براساس آن ضریب پواسون نسبت‌به دو 


سال سی و سوم شمارة سهء ۱۳۹۹ 


۳۷ تحلیل احتمالاتی نشس تآنی شالوده‌ی مستطیلی بااستفاده از نرمافزارهای... 


9۹ 


5 
۳ 7 
4465 می‌مترداا م 
) 2 
با ۴۷۵۱0330 ۳۵۲ 
0 ید۷ 


شکل ۸ نمودار توزیع فراوانی احتمال شکست در نرم‌افزار 685 


بردار حساسیت آلفا بردار حساسیت اتا بردار حساسیت دلتا 


شکل ٩‏ نمودار بردارهای حساسیت متغیرهای تصادفی 


شکل ۱۰ نمودار بردارهای حساسیت متغیرهای تصادفی 


نشریه مهندسی عمران فردوسی سال سی و سوم شماره سه. ۱۳۹۹ 


۳/۸ 


تحلیل احتمالاتی نشست آنی شالوده مستطیلی بااستفاده از نرمافزارهای... 


جدول ۶ بردارهای حساسیت متغیرهای تصادفی 


نوع پردار حساسیت بار وارده 
۳ ۰۸۱/۳۷۸ 
7 ۲ 
0 ۲ 
8 ۷۸ 


۳۱/۰۸ 5 


0 م0۷ ۵0 ۶۳۲2 5۷ 8۵0۷۵۵ عنا108 


ون 1۳2۲ 5 


5 ۷۰ 


25 


ضریب پواسون مدول الاستیسیته 
۱ ۲۸/- 
- ۰ - 
۸/۲ ۰/۱۱۷۱ 
,"+ 2 
۲- ۲- 
9۹ 


0۹ 


شکل ۱۱ تحلیل حساسیت متغیرهای تصادفی برمبنای قابلیت اعتماد 


نتیجه گیری 
در این تحقیق تحلیل قابلیت اعتماد نشست آنی یک 
شالوده‌ی مستطیلی بااستفاده از روش‌های مرتبه‌اول و 
مرتبه‌دوم قابلیت اعتماد و شبیه‌سازی مونت‌کارلو انجام 
شد. بررسی‌ها نشان داد که تحلیل‌های انجام‌شده با 
روش‌های نام‌برده با اختلاف نتایج بسیار کمی, بالا بودن 
انختمال شکست تین اه مرا تسیک شالو خی 
سطحی موردنظر را نسبت به احتمال شکست قابل قبول 
در زمینه‌ی مهندسیء تأیید می‌کند. هم‌چنین مقایسه‌ی 
احتمال شکست محاسبه‌شده از روش‌ها نشان داد باتوجه 
به نتیجه‌ی مشابه روش مرتبه‌دوم قابلیت اعتماد به 
شبیه‌سازی مونت‌کارلو. که دقت بسیار بالایی در تخمین 
احتمالات دارد. رویه‌ی حالت حدی (0 < (0)6) به 
سهمی نزدیک‌تر است. به همین دلیل روش مرتبه‌اول با 
تقریب خطی از رویه‌ی حالت حدی را با خطا مواجه 


نشریه مهندسی عمران فردوسی 


می‌کند. به این علت که در شبیه‌سازی مونت کارلو برای 
کاهش تغییرات در تخمین احتمال شکست و رسیدن به 
جواب با دقت بالا نیاز به تولید تعداد زیادی داده‌ی 
تصادفی و تکرار شبیه‌سازی (280000 < 0) وجود دارد. 
برای مسئله‌ی حاضر روش مرتبه‌دوم قابلیت اعتماد با 
هزینه‌ی محاسباتی بسیار کمتر. جایگزین مناسبی برای 
شبیه‌سازی مونت کارلو می‌باشد. 

تحلیل حساسیتی که درمورد متغیرهای تصادفی و 
بااستفاده از بردارهای مختلف انجام شد. متغیرها را ازنظر 
آشبنه .۵ ال گلاری: نر. یزان تفس شالوهه 
اولویت‌بندی کرده‌است. بار واردشده بر شالوده‌ی سطحی 
دارای بیشترین اهمیت است و تغییر آن تأثبر زیادی بر 
نشست دارد. پس از این متغیر ضریب پواسون که از 
خصویات محیطی مسئله به‌شمار می رود دارای اهمیت 


سال سی و سوم شمارة سهء ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی- سام مهرنیا 


محیط مورد نظر احتمالات محاسبه‌شده در باب ایمنی را پهبود می‌بخشد. 
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